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During  a  single-string  integration  test  of  the  NASA  HiVHAc  Hall  thruster,  a  number  of 
plasma  diagnostics  were  implemented  to  study  the  effect  of  varying  facility  background 
pressure  on  thruster  operation.  These  diagnostics  characterized  the  thruster  performance, 
the  plume,  and  the  plasma  oscillations  in  the  thruster.  Thruster  performance  and  plume 
characteristics  as  functions  of  background  pressure  were  previously  published.  This  paper 
will  focus  on  changes  in  the  plasma  oscillation  characteristics  with  changing  background 
pressure.  The  diagnostics  used  to  study  plasma  oscillations  include  a  high-speed  camera  and 
a  set  of  high-speed  Langmuir  probes.  The  results  show  a  rise  in  the  oscillation  frequency  of 
the  ‘‘breathing”  mode  with  rising  background  pressure,  which  is  hypothesized  to  be  due  to  a 
shortening  acceleration/ionization  zone.  An  attempt  is  made  to  apply  a  simplified  ingestion 
model  to  the  data.  The  combined  results  are  used  to  estimate  the  maximum  acceptable 
background  pressure  for  performance  and  wear  testing. 


Abbreviations 

GRC  =  Glenn  Research  Center  MCD 

SMD  =  Science  Mission  Directorate  CEX 

ISPT  =  In-Space  Propulsion  Technology  SEE 

HiVHAc  =  High  Voltage  Hall  Accelerator  Cl 

EDU  =  Engineering  Development  Unit 


I.  Introduction 

HIS  paper  is  a  continuation  of  a  previous  work  on  the  effect  of  background  pressure  on  the  performance  and 
plume  of  the  NASA  High-Voltage  Hall  Accelerator  (HiVHAc)  thruster/  While  the  previous  paper  focused  on 
the  breakdown  of  performance  into  various  efficiency  factors  and  far-field  time-averaged  diagnostics,  this  paper  will 
focus  on  plasma  oscillation  characteristics  and  time -resolved  diagnostics. 

The  HiVHAc  project  is  funded  by  the  NASA  Science  Mission  Directorate’s  (SMD’s)  In-Space  Propulsion 
Technology  (ISPT)  program  with  the  goal  of  delivering  a  Hall  thruster  system  beyond  the  current  state-of-the-art 
(SotA).  The  project  objective  is  to  enable  various  Discovery  class  missions  while  keeping  the  cost  of  the  propulsion 
system  low.  Mission  studies  in  2009  showed  that  a  3.5-kW  HiVHAc  is  able  to  outperform  an  SotA  4.5-kW  flight 
Hall  thruster  by  delivering  6-12%  more  payload  mass  on  four  missions  of  interest.^  These  missions  include  Vesta- 
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=  Confidence  Interval 
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Ceres  rendezvous  (the  Dawn  mission),  Koppf  comet  rendezvous,  Nereus  (a  near-Earth  asteroid)  sample  return,  and 
NEARER,  which  involves  two  near-Earth  asteroid  returns.  The  thrust  profiles  of  the  missions  range  from  ones  that 
favor  high  specific  impulse  to  ones  that  favor  high  thrust-to-power  demonstrating  the  wide  throttle-ability  of  the  3.5- 
kW  HiVHAc. 

The  iteration  of  the  HiVHAc  Hall  thruster  tested  in  this  work  is  the  Engineering  Development  Unit  2  (EDU2), 
which  is  designed  to  operate  at  up  to  3.9  kW  discharge  power.  This  version  of  the  HiVHAc  features  a  discharge 
channel  replacement  mechanism  designed  to  guarantee  long  operational  life  time  for  high  specific  impulse  (high 
discharge  voltage)  operation.  Compared  to  the  previous  iteration,  which  is  the  EDU 1 ,  the  EDU2  has  better  thermal 
management,  simpler  design  for  the  channel  replacement  mechanism,  and  superior  voltage  isolation.  The  EDU2  has 
undergone  the  performance  acceptance  test  (PAT),^  a  vibration  test,^  a  plume  divergence  characterization  test,"^  an 
in-situ  mechanism  actuation  tests  (unpublished),  and  a  single-string  integration  test  (SSIT).^  The  background 
pressure  study  described  in  this  paper  was  performed  during  the  SSIT  in  April  and  May  of  2013.^  This  work 
produced  data  that  can  guide  future  HiVHAc  testing,  including  the  wear  test. 

To  reach  an  environment  as  close  to  space-like  as  possible,  the  EDU2  was  tested  in  NASA  Glenn  Research 
Center’s  (GRC’s)  Vacuum  Facility  5  (VF5).  Many  diagnostics  were  deployed  as  a  part  of  the  testing,  including  a 
thrust  stand,  a  retarding  potential  analyzer,  a  Wien  filter,  also  known  as  an  ExB  probe,  a  Faraday  probe,  an  infrared 
camera,  a  high-speed  camera,  and  a  set  of  high-speed  Langmuir  probes.  Additional  details  regarding  the  SSIT  test 
setup  can  be  found  in  Kamhawi’s  paper. ^  The  present  paper  will  focus  on  setup  and  analysis  of  the  high-speed 
camera  and  the  high-speed  Langmuir  probes,  as  well  as  the  implications  of  the  background  pressure  study  for  future 
testing. 


II.  Past  Work 

Many  kinds  of  plasma  oscillations  have  been  known  to  exist  in  Hall  thrusters.^  Since  the  1960’s,  many  studies 
have  been  performed  to  characterize  these  oscillations.^'^®  Yet  the  roles  of  these  oscillations  in  relation  to  thruster 
performance  and  life  time  are  still  largely  unknown.  Currently,  there  are  two  oscillation  modes  that  are  believed  to 
be  important  to  Hall  thruster  operations,  the  “breathing”  mode  and  the  “spokes”  mode.  Both  modes  are  believed  to 
be  tied  to  ionization  processes  in  the  Hall  thruster  because  their  frequencies  are  in  the  kilohertz  to  tens  of  kilohertz 
range,  which  is  most  closely  related  to  the  characteristic  ion  and  neutral  transit  time  across  the  discharge  axial  length 
scale.®’  Recently,  there  is  a  new  surge  of  interest  in  these  plasma  oscillations  because  they  have  been  observed  to 
vary  greatly  with  facility  conditions'^  and  may  drive  differences  between  thruster  behavior  in  ground  tests  and  on 
orbit.  As  more  Hall  thrusters  are  being  considered  for  flight  applications,  including  the  HiVHAc  series,  it  has 
become  critical  for  the  Hall  thruster  community  to  better  understand  the  behavior  of  these  plasma  oscillations  so  that 
on-orbit  performance  can  be  predicted  using  ground  test  data.  For  the  present  study,  the  high-speed  camera  and 
Langmuir  probe  were  chosen  to  study  said  plasma  oscillations. 

The  high-speed  camera  has  only  recently  been  developed  into  a  quantitative  tool  for  studying  Hall  thruster 
plasma  oscillations.  High-speed  optical  imaging  on  a  Hall  thruster  was  first  performed  by  Darnon,  et.  al.,  though  the 
view  was  not  perpendicular  to  the  annulus  plane  and  did  not  capture  azimuthal  plasma  movement.  The  instrument 
Darnon,  et.  al.,  used  was  a  charge-coupled  device  with  an  image  intensifier.^^  Subsequently,  Lobbia,  Liu,  and 
Gallimore  studied  a  cluster  of  annular  Hall  thrusters  using  a  commercial  high-speed  camera^^  while  Parker,  Raitses, 
and  Fisch  studied  a  cylindrical  Hall  thruster.^^  Shortly  after,  McDonald  developed  a  Fourier-transform  technique  to 
study  both  breathing  and  spokes  mode  in  annular  Hall  thrusters.^®  Sekerak  then  combined  the  Fourier-transform 
technique  with  results  from  high-speed  Langmuir  probes  to  study  oscillation  mode  transition  in  Hall  thrusters  as  a 
function  of  magnetic  field  strength.^®  The  analysis  technique  used  in  the  present  study  is  based  on  the 
aforementioned  Fourier-transform  technique. 

The  high-speed  Langmuir  probe  was  first  used  on  a  Hall  thruster  by  Lobbia  and  Gallimore,  who  first  identified 
all  of  the  distortion  effects  associated  with  high-speed  sweeping  of  Langmuir  probes  up  to  the  megahertz  range.^^ 
Lobbia  then  combined  a  high-speed  dual  Langmuir  probe  (HDLP)  with  a  Fourier-transform  correlation  technique  to 
temporally  and  spatially  resolve  the  breathing  mode  in  a  Hall  thruster.  The  high-speed  Langmuir  probe  technique 
used  by  Sekerak  is  a  refined  version  of  Lobbia’s  work.^® 

The  breathing  mode  oscillation  is  believed  to  be  an  axial  ionization  wave.  The  name  comes  from  the 
resemblance  of  the  exhausted  plasma  packets  to  that  of  breaths  of  gas  during  breathing  mode  oscillation.  The  current 
best  quantitative  description  of  the  breathing  mode  oscillation  is  based  on  modeling  work  done  by  Fife.  Fife  started 
with  the  assumption  that  the  breathing  mode  is  an  axial  ionization-driven  oscillation  and  derived  the  scaling 
relationship  shown  in  Eq.  (1)  below. 
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In  Eq.  (1),  fi  is  the  frequency  of  the  breathing  mode,  Ui  is  the  speed  of  the  ions  leaving  the  thruster,  Un  is  the  speed  of 
the  neutrals  entering  the  ionization  zone,  and  Li  is  the  characteristic  length  over  which  the  oscillation  action  is  taking 
place.  Typically,  Li  is  taken  to  be  the  overall  length  of  the  discharge  zone,  where  both  ionization  and  acceleration  are 
taking  place.  This  scaling  relationship  has  been  shown  to  predict  general  trends  in  the  breathing  mode  oscillation 
frequency  with  discharge  voltage  though  the  exact  values  do  not  always  agree.  The  role  of  the  breathing  mode  in 
electron  transport  is  currently  unclear,  though  Lobbia  has  argued  that  the  breathing  mode  is  an  intrinsic  feature  in  the 
Hall  thruster  and  is  necessary  for  the  thruster  to  achieve  enough  electron  transport  to  maintain  the  discharge.  The 
scaling  relation  was  not  created  with  the  ingestion  of  background  neutrals  in  mind. 

The  spokes  mode  oscillation  is  believed  to  be  an  azimuthal  ionization  wave.  The  name  comes  from  the 
appearance  of  rotating  plasma  packets,  which  resemble  spokes  on  a  wheel,  in  the  channel  of  a  Hall  thruster  when 
viewed  face-on.  The  current  best  understanding  of  the  spokes  mode  oscillation  come  from  works  by  Janes  and 
Lowder^  and  McDonald.  The  spokes  mode  is  presented  as  an  alternate  means  by  which  electron  current  can  be 
transported  across  magnetic  field  lines.  Subsequent  work  by  Sekerak  showed  that  a  Hall  thruster  can  operate  with 
almost  no  breathing  mode  and  strong  spokes  mode  oscillations  present.  Sekerak  observed  that  the  presence  of 
spokes  mode  without  breathing  mode  is  correlated  with  high  thrust  efficiency  for  a  medium-power  Hall  thruster. 
Spokes  have  been  observed  to  travel  in  the  same  direction  as  the  Hall  current  and  anywhere  from  two  to  six  spokes 
can  be  observed  simultaneously.^^’  There  are  currently  no  existing  scaling  relationships  for  the  number  and  speed 
of  the  spokes  present. 


III.  Experimental  Setup 

The  experimental  setup  was  previously  reported^  and  is  repeated  here  with  new  details  about  the  high-speed 
diagnostics  for  convenience.  To  simplify  plot  labeling,  operating  conditions  will  be  labeled  as  vvv-k.k-Pnx,  where 
vvv  is  the  discharge  voltage  in  volts,  k.k  is  the  discharge  power  in  kilowatts,  and  n  is  the  normalized  background 
pressure,  defined  as  the  number  of  multiples  of  the  lowest  achievable  operating  background  pressure  for  the  throttle 
point  vvv-k.k.  This  pressure  is  measured  by  an  ion  gauge  located  right  below  the  thrust  stand.  Lor  example  Plx 
indicates  the  operating  condition  corresponding  to  the  minimum  achieved  background  pressure,  while  PI  Ox 
indicates  the  background  pressure  is  ten  times  that  of  Plx.  Unless  otherwise  noted,  all  spatial  positions  presented  in 
this  paper  have  been  normalized  by  the  mean  channel  diameter  (MCD)  of  the  thruster.  MCD  is  defined  as  the 
average  of  the  inner  and  outer  discharge -channel  wall  diameters. 

A.  Test  Article  and  Test  Matrix 

The  NASA  HiVHAc  EDU2  is  a  3.9-kW  xenon  Hall 
thruster.  The  thruster  is  highly  throttle-able.  The  discharge 
voltage  can  vary  from  200  to  650  V,  corresponding  to  specific 
impulse  values  of  1200  to  2700  s.  The  discharge  power  ranges 
from  0.3  to  3.9  kW.  During  testing,  the  magnetic  field  settings 
are  selected  to  create  a  magnetic  lens  symmetric  about  the 
channel  centerline,  while  maximizing  total  efficiency.  This 
magnetic  field  setting  optimization  is  performed  at  the  lowest 
achievable  background  pressure  for  each  throttle  point.  Ligure  1 
shows  a  picture  of  the  NASA  HiVHAc  EDU2  with  various  test 
equipment. 

Lor  the  purpose  of  this  test,  seven  throttling  points  were 
selected.  These  seven  points  are  300-1.5,  300-3.0,  400-2.0,  400- 
3.2,  500-2.0,  500-3.0,  and  500-3.9.  At  each  throttling  point,  the 
thruster  settings  (including  magnet  currents)  are  optimized  at 
the  lowest  achievable  pressure  as  measured  by  the  ion  gauges 
next  to  the  thrust  stand.  Then,  the  background  pressure  was 
raised  by  injecting  xenon  via  an  auxiliary  flow  line  that  exits  at 
>4  meters  downstream  of  the  thruster  pointed  away  from  the 
thruster. 


Figure  1.  NASA  HiVHAc  EDU2  and  thrust 
stand  setup. 
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B.  Test  Facility 

Testing  was  conducted  in  NASA  GRC’s  VF5.  This  cylindrical  facility  is  4.6  m  in  diameter,  18.3  m  long,  and  is 
pumped  with  a  set  of  cryo-panels  and  20  oil  diffusion  pumps.  The  thruster  was  mounted  on  a  thrust  stand  located 
close  to  the  cryo-panels,  with  the  thruster  firing  away  from  the  panels.  Facility  pressures  were  monitored  with  four 
ion  gauges,  three  of  which  were  mounted  next  to  the  thrust  stand  and  the  fourth  mounted  on  the  facility  wall  at  the 
same  axial  location  as  the  thruster.  Manufacturer  specifications  state  that  the  ion  gauges  are  accurate  to  ±6%  of 
reading.  The  positions  of  the  gauges  are  shown  in  Fig.  1.  Ion  gauges  1  and  2  are  both  facing  downstream  while  ion 
gauge  3  is  facing  upstream.  Ion  gauge  1  and  2  agree  to  within  10%  of  each  other.  Ion  gauge  3  reports  0.63  to  0.72 
times  the  reading  of  ion  gauge  2.  Ion  gauge  2  readings  were  used  to  determine  the  number  of  multiples  of  the  lowest 
achievable  background  pressure  that  the  thruster  was  experiencing.  The  lowest  pressure  achieved  for  each  tested 
condition,  corrected  for  xenon,  varied  from  1.4x10'^  to  2.2x10'^  Torr. 

Thruster  propellant  was  supplied  via  commercially  available  mass  flow  controllers.  These  mass  flow  controllers 
were  calibrated  using  xenon  prior  to  testing.  Typical  uncertainty  of  measurement  was  ±1.0%  of  reading.  Research- 
grade  xenon  was  used. 

Thruster  electrical  power  was  supplied  via  the  Colorado  Power  Electronics  (CPE)  brassboard  #2  (BB#2)  power 
processing  unit  (PPU).^"^  This  PPU  is  capable  of  supplying  200  to  700  V  at  up  to  4  kW  over  the  main  discharge  line 
as  well  as  supply  power  to  the  cathode  heater,  keeper,  and  magnet  electrical  circuits.  The  PPU  was  located  outside 
of  the  vacuum  facility  during  testing. 


C.  High-Speed  Diagnostics 

The  high-speed  camera  was  a  commercial  device  on  loan  from  the  Aerospace  Corporation.  To  obtain  a  face-on 
view  of  the  thruster  without  putting  the  high-speed  camera  into  the  plume  a  large  plane  mirror  was  deployed.  The 
camera  was  positioned  outside  of  the  vacuum  chamber  and  views  the  mirror  through  a  window  port.  The  camera 
was  set  to  capture  at  a  frame  rate  of  roughly  80,000  frames  per  second  at  256  by  256  pixels.  In  order  to  capture  the 
cathode  as  well  as  compensate  for  possible  misalignment  during  pump  down  the  view  area  is  quite  a  bit  larger  than 
the  area  taken  up  by  the  discharge  channel,  which  took  up  a  100-by- 100-pixels  area.  Eigure  2  shows  an  overhead 
diagram  of  the  positions  of  the  mirror  and  camera.  The  view  of  the  camera  is  actually  4°  off  to  the  side  instead  of 
perfectly  straight-on  due  to  the  need  to  accommodate  an  IR  camera.  This  small  offset  was  considered  acceptable  as 
it  produces  negligible  impact  on  the  data  analysis.  Due  to  computer  error,  high-speed  video  at  varying  background 
pressure  is  only  available  for  the  throttle  points  500-2.0,  500-3.0,  and  500-3.9.  High-speed  video  was  obtained  at 
only  the  lowest  achievable  pressure  for  the  throttle  points  300-1.5,  400-2.0,  and  400-3.2. 

The  high-speed  Langmuir  probe  array  consists  of  four  pairs  of  Langmuir  probes.  Each  pair  consists  of  one  active 
probe  and  one  null  probe.  Data  from  the  null  probe  is  used  to  correct  for  capacitance  effects,  which  allows  for 
accurate  high-speed  Langmuir  probe  sweeps.  The  entire  probe  array  is  mounted  on  a  two-axis  motion  system  and  is 
well  outside  the  plume  when  not  in  use.  The  four  pairs  of  probes  are  divided  into  two  groups.  The  upstream  group 
consists  of  two  pairs  set  2  cm  apart  in  the  azimuthal  direction.  The  downstream  group  consists  of  the  other  two  pairs 
set  6  cm  apart  in  the  azimuthal  direction.  The  axial  distance  between  the  upstream  and  downstream  groups  is  10  cm. 
There  are  three  measurement  positions  for  the  array  at  each  thruster  condition.  Eor  these  three  positions,  the 
upstream  probe  is  10,  15,  and  20  cm  away  from  the  thruster  exit  plane,  respectively,  and  the  probe  tips  are  radially 
aligned  with  the  centerline  of  the  channel.  To  correlate  the  Langmuir  probe  signals  across  different  locations,  the 
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Figure  2.  A  diagram  of  the  diagnostics  setup. 
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Figure  3.  A  diagram  of  the  movement 
path  of  the  high-speed  Langmuir 
probe  array. 
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discharge  electrical  line  is  monitored.  Figure  3  shows  an  overhead  diagram  of  the  movement  path  of  the  probe  array 
and  the  interrogation  locations  with  respect  to  the  thruster. 

A  specialized  high  frequency  power  supply  is  used  to 
supply  a  200-kHz  signal  through  custom  circuitry  to  sweep 
the  probes.  The  custom  circuitry  is  tuned  to  minimize 
capacitive  effects  associated  with  high-speed  probe 
sweeping.  A  custom  computer  with  high-speed  data 
acquisition  systems,  high  processing  power,  and  high 
storage  capacity  is  used  to  record  the  data.  Figure  4  shows  a 
photograph  of  the  high-speed  Langmuir  probe  array  moving 
into  place  to  record  data  while  the  thruster  was  operating. 

IV.  Data  Analysis 

A.  High-Speed  Camera  Analysis 

For  this  study,  the  data  analysis  focused  on  the  breathing 
and  spokes  mode  oscillations.  The  high-speed  camera 
videos  are  processed  in  the  following  steps. 

1 .  The  image  is  cropped  so  that  only  the  annulus  region  is 
analyzed. 

2.  The  video  is  averaged  over  time  to  create  an  averaged  image. 

3.  A  circle  fit  is  performed  on  the  averaged  image  to  find  the  center  of  the  thruster  and  the  boundaries  of  the 
discharge. 

4.  The  data  is  divided  into  120  azimuthal  bins  and  the  average  intensity  in  each  bin  on  each  frame  is  calculated. 

5.  The  binned  data  is  normalized  against  the  time-averaged  values  of  the  binned  data. 

6.  A  two-dimensional  (2D)  Fourier  transform  is  performed  to  obtain  a  series  of  power  spectra  categorized  by  the 
mode  number  m. 

The  data  resulting  from  step  5  are  basically  the  relative  fluctuation  of  each  azimuthal  bin  with  time.  This  step 
produces  results  that  are  much  less  sensitive  to  variations  in  pixel  sensitivity  to  light,  transmission  of  the  viewport, 
and  reflectance  of  the  mirror.  Any  breathing  or  spokes  mode  present  will  show  up  more  clearly  when  this  step  is 
performed.  Note  that  step  5  also  removes  any  real,  time-averaged,  azimuthal  variation  in  light  intensity.  For 
instance,  the  presence  of  a  stationary  hot  spot  would  be  removed  by  this  step.  This  shortfall  is  considered  acceptable 
because  the  present  study  is  focused  on  oscillatory  instead  of  stationary  features.  For  the  power  spectra  obtained 
from  step  6,  the  m  =  0  spectrum  shows  the  frequency  of  any  breathing  mode  present,  the  m  =  1  spectrum  shows  the 
frequency  of  any  spokes  mode  with  one  spoke  present,  the  m  =  2  spectrum  shows  the  frequency  of  any  spokes  mode 
with  two  spokes  present,  and  so  on. 

To  properly  interpret  the  power  spectra  from  the  2D 
Fourier  transform,  one  must  first  understand  what  the 
transform  is  really  doing  to  the  data.  The  transform  breaks 
down  the  data  into  a  series  of  2D  sinusoids  that  when  taken 
together  can  re-constitute  the  original  data.  Since  one  of  the 
dimensions  of  the  binned  data  is  time  and  the  other  is 
azimuthal  angle,  the  transform  looks  for  waves  that  are  both 
repeating  in  time  and  in  angular  space.  Figure  5  shows 
illustrations  of  the  ideal  modal  patterns  that  the  transform  is 
looking  for  when  processing  the  data.  Parts  of  the  annulus 
with  higher  than  average  light  intensity  are  shown  in  white 
while  parts  with  lower  than  average  light  intensity  are 
shown  in  black.  Although  the  true  ideal  pattern  will  have  an 
intensity  that  varies  in  a  sinusoidal  fashion  as  opposed  to 
having  discrete  black  and  white  shades,  this  simplified 
illustration  is  used  for  better  clarity. 

An  important  feature  of  the  modal  patterns  in  Fig. 5  is  that  the  angular  length  of  the  light  packet  is  necessarily 
shorter  for  higher  mode  number  patterns  than  for  lower  mode  number  patterns.  For  example,  the  light  packets  of  m 
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B'igure  5.  Illustrations  of  ideal  modal  patterns  as 
ipplied  in  the  2D  Fourier  transform  analysis. 


Figure  4.  A  photograph  of  the  high-speed 
Langmuir  probe  array  with  the  thruster 
operating.  Only  the  downstream  group  of  probes 
is  visible  due  to  the  camera  angle. 


=  3  mode  are  60  degrees  long  while  the  light  packet  of  m  =  1  mode  is  180  degree  long.  Another  important  feature  is 
that  these  model  patterns  have  light  packets  that  are  equally  spaced.  However,  real  Hall  thruster  plasma  waves  are 
generally  more  chaotic,  with  light  packets  that  can  be  shorter  or  longer  than  ideal  and  with  spacing  that  are  not  even. 
Thus,  a  spokes  mode  oscillation  pattern  with  two  short  light  packets  may  be  de-composed  into  a  series  of  modal 
patterns  at  the  m  =  4  mode  in  addition  to  the  m  =  2  mode  because  the  length  of  the  light  packets  are  closer  to  that  of 
the  ideal  m  =  4  modal  pattern  but  the  spacing  of  the  packets  are  more  similar  to  the  ideal  m  =  2  modal  pattern. 
Another  example  is  when  4  light  spokes  are  present  with  uneven  lengths  and  spacing,  and  is  decomposed  into  a 
series  of  peaks  at  m  =  3,  4,  and  5.  The  real  spokes  pattern  did  not  fit  neatly  into  the  ideal  modal  patterns  but  that 
does  not  necessarily  mean  there  were  3,  4,  and  5  light  spokes  present  at  the  same  time;  instead,  the  real  pattern  is  a 
super-positioning  of  m  =  3,  4,  and  5  modal  patterns  on  top  of  each  other.  To  complicate  matters  further,  if  the 
thruster  changes  from  an  ideal  m  =  3  modal  pattern  to  an  ideal  m  =  4  modal  pattern,  and  then  changes  to  an  ideal  m 
=  5  modal  pattern  during  the  filming  of  the  high-speed  video,  the  resulting  Fourier  transform  can  also  produce  a 
series  of  peaks  at  m  =  3,  4,  and  5.  Thus,  close  examination  of  the  actual  high-speed  video  is  sometimes  necessary  to 
determine  the  exact  nature  of  the  modal  pattern  present. 

In  contrast,  analysis  of  the  breathing  mode  oscillation  can  be  completed  by  either  averaging  the  light  intensity 
data  across  all  bins  and  performing  a  ID  Fourier  transform  or  taking  the  m  =  0  result  from  the  2D  Fourier  transform. 
Both  approaches  are  performed  in  the  present  analysis  to  check  for  consistency. 

B.  High-Speed  Langmuir  Probe  Analysis 

High-speed  Langmuir  probe  analysis  is  described  in  much  greater  detail  in  Lobbia’s  work^^  and  will  only  be 
briefly  summarized  here.  The  high-speed  Langmuir  probe  data  is  first  corrected  for  capacitive  effects  associated 
with  high-speed  sweeping.  This  is  done  by  subtracting  the  null  probe  signal  from  the  active  probe  signal.  The  data  is 
then  chopped  into  individual  sweeps  and  analyzed  using  standard  collisionless  non-magnetized  thin-sheath 
Langmuir  probe  theory.  Fourier  transforms  are  performed  on  the  resulting  number  density  as  a  function  of  time  to 
obtain  power  spectra. 


V.  Results 


A.  High-Speed  Camera  Results 

The  discharge  of  the  HiVHAc  EDU2  was  found  to  be  dominated  by  breathing  mode  oscillations.  For  the  500  V 
throttle  points  at  all  test  background  pressures,  there  were  no  discernable  signs  of  spokes  mode  oscillations.  Spokes 
were  found  at  other  throttle  points  but  were  found  to  exist  only  within  limited  sectors  around  the  thruster. 

Figures  6-9  show  the  theta-t  diagrams  (also  called  spokes  surfaces)  for  the  500-3. 9-Pnx  operating  conditions. 
The  theta-t  diagram  is  a  surface  contour  plot  where  the  x-axis  is  the  azimuthal  axis  and  the  y-axis  is  the  time  axis. 
Red  indicates  higher-than-average  light  intensity  and  blue  indicates  lower-than-average  light  intensity.  Theta,  the 
azimuthal  coordinate,  is  defined  as  0  degrees  at  the  top  of  the  annulus.  Increasing  theta  is  equivalent  to  clockwise 
movement  on  the  image.  Note  the  video  is  filmed  through  a  mirror  so  the  filmed  image  is  a  horizontally-mirrored 
image  of  the  real  thruster  face  and  the  reported  theta  coordinate  is  not  corrected  for  this  reflection. 
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Figure  6.  Theta-t  diagram  for  500-3.9-Plx. 
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Figure  7.  Theta-t  diagram  for  500-3.9-P3x. 
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Figure  8.  Theta-t  diagram  for  500-3.9-P5x. 
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Figure  9.  Theta-t  diagram  for  500-3.9-P10x. 


From  Figs.  6-9,  one  can  see  a  large  number  of  horizontal  bands  of  color.  These  bands  of  color  indicate  a  fairly 
uniform  light  intensity  across  all  azimuthal  locations  that  varies  in  time,  which  is  the  signature  of  breathing  mode 
oscillations.  In  contrast,  a  spokes  oscillation  shows  up  as  bands  of  diagonal  lines  because  the  light  packet  is  traveling 
across  both  angle  and  time.  Examination  of  500-2.0  and  500-3.0  throttle  point  data  as  functions  of  background 
pressure  indicates  the  same  trend;  there  is  a  strong  breathing  mode  of  varying  frequency  present  but  no  spokes  mode 
oscillations. 

Figures  10-13  show  the  high-speed  camera  power  spectra  for  the  500-3. 9-Pnx  operating  conditions  with  mode 
number  0  to  5  being  shown.  Note  m  =  0  corresponds  to  breathing  mode  while  m  =  l-i-  correspond  to  spokes  mode 
with  different  modal  patterns.  The  y-axis  of  the  figures  is  spectral  power  on  a  logarithmic  scale  and  x-axis  is 
frequency  in  kilohertz.  Since  the  videos  were  taken  at  80,000  frames  per  second,  the  maximum  resolved  wave 
frequency  is  40  kHz. 


Frequency,  kHz 

Figure  10.  High-speed  camera 
power  spectra  for  500-3.9-Plx. 
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Figure  11.  High-speed  camera 
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Figure  12.  High-speed  camera 
power  spectra  for  500-3.9-P5x. 


Figure  13.  High-speed  camera 
power  spectra  for  500-3.9-P10x. 


While  some  peak  structures  can  be  seen  for  m  =  1  and  2  the  peak  magnitude  of  these  structures  are  two  orders  of 
magnitude  lower  than  the  m  =  0  peak  structures.  The  breathing  mode  oscillation  is  clearly  dominant  over  spokes 
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mode  oscillation  at  these  operating  conditions.  Analysis  of  the  500-2.0  and  500-3.0  operating  conditions  show  the 
same  trend. 

Although  high-speed  videos  for  the  background  pressure  test  at  the  300  V  and  400  V  throttle  points  are  not 
available,  there  are  videos  at  the  lowest  achievable  pressure  that  were  obtained  during  the  performance  testing  phase 
of  the  single-string-integration  test.  A  search  through  these  videos  was  performed  to  help  determine  why  spokes 
mode  oscillation  was  absent  in  the  videos  of  the  500  V  operation.  Examination  of  thruster  operation  at  400-3. 2-P lx 
revealed  a  peculiar  phenomenon  where  spokes  mode  is  observed  but  only  over  a  small  section  of  the  thruster 
annulus.  Figures  14  and  15  show  the  theta-t  diagrams  and  the  high-speed  camera  power  spectra,  respectively,  for  the 
400-3. 2-P  lx  operating  condition. 


Theta,  deg. 

Figure  14.  Theta-t  diagram 
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Figure  15.  High-speed  camera  power  spectra 
for  400-3.2-Plx. 


Figure  14  indicates  that  a  single  spoke  repeatedly  appears  and  disappears  in  the  space  between  theta  equal  to  90 
and  135  degrees.  The  camera  power  spectra  appear  to  indicate  that  the  size  of  the  spoke  is  a  good  match  for  the  m  = 
4  and  5  modal  patterns.  The  frequency  of  the  spokes  is  distinctly  different  from  that  of  the  breathing  mode 
oscillation.  Direct  observation  of  the  high-speed  video  confirms  the  presence  of  this  peculiar  spokes  mode  structure. 
Spokes  mode  structure  does  appear  to  form  within  the  HiVHAc  EDU2  but  some  damping  mechanism  damps  out 
these  spokes  shortly  after  formation.  Examination  of  the  corresponding  sector  of  the  thruster  post-test  did  not  reveal 
any  obvious  differences  compared  to  the  rest  of  the  thruster.  While  it  is  still  unclear  why  spokes  mode  oscillations 
could  not  be  sustained  in  the  HiVHAc  EDU2  operating  under  the  tested  conditions,  it  is  clear  that  the  analysis  is 
capable  of  capturing  spokes,  and  that  spokes  oscillations  can  be  formed.  The  spokes  oscillations  are  just  not  being 
sustained. 

Table  1  summarizes  the  breathing  mode  frequency  and  magnitude  as  functions  of  facility  background  pressure 
calculated  from  the  high-speed  videos.  Oscillation  magnitude  is  reported  in  the  form  of  relative  deviation  in  light 
intensity  with  a  95% -confidence-interval  (95%-CI),  much  like  how  relative  uncertainty  is  typically  reported.  The 
relative  deviation  is  reported  as  a  percentage  of  the  time-averaged  light  intensity,  where  95%  of  the  time  the  light 
intensity  falls  within  the  reported  boundary.  For  example,  if  50%  relative  deviation  is  reported,  that  means  95%  of 
the  time,  the  light  intensity  fell  within  ±50%  of  the  average  light  intensity.  This  table  shows  that  the  breathing  mode 
frequency  rises  while  the  relative  deviation  of  light  intensity  falls  with  rising  background  pressure.  Recall  from  Eq. 
(1)  the  breathing  mode  frequency  is  theorized  to  scale  with  the  square -root  of  the  ion  and  neutral  velocity  and  the 

T able  1.  Summary  of  breathing  mode  frequency  calculated  from  high-speed  imaging  data. 


Disch. 

voltage, 

V 

Disch. 

power, 

kW 

Breathing  mode  frequency, 
kHz 

Relative  deviation  of  light 
intensity,  95%-CI 

lx 

3x 

5x 

lOx 

lx 

3x 

5x 

lOx 

300 

1.5 

9.4 

- 

- 

- 

16% 

- 

- 

- 

400 

2.0 

13.0 

- 

- 

- 

13% 

- 

- 

- 

400 

3.2 

20.4 

- 

- 

- 

14% 

- 

- 

- 

500 

2.0 

22.0 

24.0 

25.7 

- 

47% 

28% 

25% 

- 

500 

3.0 

20.6 

29.2 

29.5 

34.4 

36% 

31% 

32% 

23% 

500 

3.9 

26.8 

28.9 

30.2 

36.2 

42% 

29% 

29% 

23% 
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inverse  of  the  characteristic  length.  At  a  given  throttle  point,  retarding  potential  analyzer  data  previously  reported 
shows  that  the  ion  velocity  remains  roughly  constant  as  background  pressure  varies.^  The  neutral  velocity  scales 
with  the  square  root  of  the  anode  temperature  and  is  expected  to  remain  roughly  constant  as  background  pressure 
varies.  The  implication  is  that  the  breathing  mode  frequency  must  be  rising  with  background  pressure  because  the 
characteristic  length  is  decreasing. 

B.  High-Speed  Langmuir  Probe  Results 

Analysis  of  the  high-speed  Langmuir  probe  data  focused  on  the  breathing  mode  oscillations.  Probe 
measurements  taken  at  different  axial  locations  are  generally  consistent  in  terms  of  measured  oscillation  frequency. 
Electron  density  and  temperature  are  observed  to  decrease  with  axial  distance  while  plasma  potential  asymptotes  to 
low  tens  of  volts.  Only  plasma  data  from  the  position  10  cm  axially  downstream  of  the  exit  plane  will  be  reported 
since  the  whole  data  set  is  too  large  to  show  here.  The  important  trends  with  background  pressure  can  be  observed 
with  only  the  10-cm-axial-position  data.  Data  from  both  of  the  probe  pairs  at  this  position  gave  the  same  answer  to 
within  measurement  uncertainty  and  the  averaged  values  are  presented. 

Table  2  summarizes  the  time-averaged  electron  density,  electron  temperature,  and  plasma  potential  at  different 
operating  conditions.  The  data  are  color-coded  by  pressure  level  for  convenience.  Generally  speaking,  the  electron 
density  increases  with  rising  background  pressure,  which  is  expected  if  ingestion  is  occurring.  Electron  temperature 
and  plasma  potential  both  decrease  with  rising  background  pressure,  suggesting  a  recession  of  the  ion  acceleration 
zone  into  the  channel  of  the  thruster. 

Table  3  summarizes  the  frequency  of  the  plasma  oscillation  and  the  magnitudes  of  oscillations  in  electron 
density  and  electron  temperature.  Oscillation  magnitude  is  reported  in  the  form  of  relative  deviation  with  a  95%- 
confidence-interval.  The  oscillation  frequencies  of  the  electron  density  and  electron  temperature  were  found  to  be 
the  same  so  only  one  frequency  is  reported.  Some  relative  deviations  of  the  electron  temperature  are  greater  than 
100%,  which  usually  happens  when  there  are  periodic  spikes  in  the  electron  temperature.  This  may  have  been  an 
artifact  of  the  data  analysis  routine  or  noise  in  the  data.  Despite  this  issue,  the  overall  qualitative  trends  in  the 
relative  deviation  of  the  electron  temperature  appear  to  be  very  consistent  across  throttle  points. 

The  oscillation  frequency  of  the  local  plasma  around  the  probe  was  found  to  be  in  fairly  good  agreement  with 
that  of  the  global  plasma  as  measured  by  the  high-speed  video.  This  agreement  gives  us  some  certainty  that  the  two 
diagnostics  are  both  interrogating  the  same  breathing  mode  phenomenon.  The  frequency  calculated  from  the  high- 


Table  2.  Summary  of  time-averaged  plasma  properties  measured  by  tbe  bigb-speed  Langmuir  probe  at  10  cm 
axially  downstream  of  tbe  exit  plane  along  tbe  tbruster  centerline. 

Discb.  Discb.  Electron  density,  m'  Electron  temperature,  Plasma  potential,  V 

voltage,  power,  _ eV _ 


V 

kW 

lx 

3x 

5x 

lOx 

lx 

3x 

5x 

lOx 

lx 

3x 

5x 

lOx 

300 

1.5 

6.4el6 

6.8el6 

- 

7.1el6 

10.1 

7.8 

- 

4.9 

53 

39 

- 

24 

300 

3.0 

1.3el7 

1.4el7 

- 

1.8el7 

5.5 

5.0 

- 

4.8 

23 

20 

- 

18 

400 

2.0 

5.3el6 

6.6el6 

- 

7.1el6 

10.6 

9.5 

- 

6.0 

54 

43 

- 

27 

400 

3.2 

8.9el6 

9.4el6 

- 

l.lel7 

7.7 

6.4 

- 

4.5 

40 

32 

- 

21 

500 

2.0 

- 

3.4el6 

3.9el6 

- 

- 

7.4 

7.0 

- 

- 

36 

32 

- 

500 

3.0 

3.0el6 

5.4el6 

6.1el6 

7.2el6 

7.7 

6.3 

5.1 

3.7 

40 

33 

27 

20 

500 

3.9 

5.4el6 

7.2el6 

8.0el6 

8.9el6 

7.5 

5.6 

4.6 

3.5 

38 

30 

25 

18 

Table  3.  Summary  of  time-resolved  plasma  properties  measured  by  tbe  bigb-speed  Langmuir  probe  at  10  cm 
axially  downstream  of  tbe  exit  plane  along  tbe  tbruster  centerline. 


Discb.  Discb.  Relative  deviation  in  Relative  deviation  in  electron  Oscillation  frequency,  kHz 
voltage,  power,  electron  density,  m'^ _ temperature,  eV _ 


V 

kW 

lx 

3x 

5x 

lOx 

lx 

3x 

5x 

lOx 

lx 

3x 

5x 

lOx 

300 

1.5 

35% 

50% 

- 

61% 

26% 

58% 

- 

127% 

10.1 

10.7 

- 

11.2 

300 

3.0 

54% 

60% 

- 

81% 

183% 

227% 

- 

303% 

29 

37 

- 

40 

400 

2.0 

37% 

44% 

- 

63% 

31% 

44% 

- 

122% 

- 

33 

- 

38 

400 

3.2 

51% 

59% 

- 

64% 

56% 

95% 

- 

212% 

30 

- 

- 

39 

500 

2.0 

- 

32% 

35% 

- 

- 

29% 

37% 

- 

- 

25 

31 

- 

500 

3.0 

59% 

34% 

37% 

41% 

23% 

36% 

47% 

138% 

22 

28 

30 

31 

500 

3.9 

38% 

37% 

43% 

53% 

27% 

40% 

72% 

220% 

30 

29 

30 

36 
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speed  video  is  of  higher  statistical  significance  because  it  was  calculated  from  data  around  the  entire  annulus  as 
opposed  to  a  single  location.  Interestingly,  the  relative  deviation  in  electron  density  and  electron  temperature 
generally  increases  with  rising  background  pressure.  This  is  the  opposite  trend  as  the  relative  deviation  in  light 
intensity  from  the  high-speed  video. 

One  possible  explanation  is  that  the  plasma  downstream  of  the  acceleration  zone  (where  the  Langmuir  probes 
are)  really  does  behave  differently  with  changing  background  pressure  than  the  plasma  inside  of  the  acceleration 
zone.  Note  that  the  high-speed  camera  signal  is  dominated  by  light  emissions  from  regions  of  the  plasma  with  high 
electron  temperature,  which  is  inside  the  ionization  and  acceleration  zones  of  the  thruster  where  electron 
temperature  reaches  tens  of  eV.  Assuming  the  behavior  of  the  plasma  inside  the  channel  differs  from  that  of  the 
plasma  outside  the  channel,  the  combined  results  of  the  high-speed  camera  and  Langmuir  probes  paint  a  very 
interesting  picture.  As  the  background  pressure  rises,  the  downstream  boundary  of  the  acceleration  zone  of  the 
thruster  is  being  pushed  upstream,  which  is  correlated  with  an  increase  in  damping  of  the  breathing  mode  magnitude 
inside  the  channel.  At  the  same  time,  breathing  mode  magnitude  outside  the  channel  becomes  less  damped.  If  we 
further  assume  that  the  breathing  mode  oscillation  enhances  electron  transport,  the  present  data  can  be  interpreted  to 
mean  that  the  shortening  of  the  ionization  and  acceleration  zones  at  high  pressure  reduces  the  amount  of  transport 
enhancement  needed  inside  the  channel  (fewer  magnetic  field  lines  to  transport  across).  At  the  same  time,  the 
increased  distance  between  cathode  and  the  downstream  boundary  of  the  acceleration  zone  at  high  pressure 
increases  the  transport  enhancement  needed  to  maintain  the  same  flow  of  electrons  to  the  channel.  While  this  theory 
is  consistent  with  the  data,  further  experiments  are  needed  to  prove  or  disprove  this  theory. 

VI.  Discussions 

A.  Correlating  Breathing  Mode  Frequency  to  Wear  Mechanics  and  Wear  Test  Recommendations 

The  present  study  is  not  able  to  form  any  conclusions  about  the  change  in  the  energy  distribution  of  the  ions 
bombarding  the  wall  of  the  discharge  channel  as  a  function  of  background  pressure.  However,  one  can  correlate  the 
breathing  mode  data  to  the  length  of  the  ionization/acceleration  zone  of  the  thruster  and  generate  recommendations 
for  testing  from  the  results. 

As  previously  stated,  assuming  the  scaling  in  Eq.  (1)  is  correct,  the  breathing  mode  frequency  only  scales  with 
ion  velocity,  neutral  velocity,  and  characteristic  length.  For  a  given  throttle  point  the  ion  and  neutral  velocities  are 
fixed  so  the  measured  breathing  mode  frequency  is  an  indirect  measure  of  the  characteristic  length.  Since  the 
characteristic  length  is  roughly  the  length  of  the  ionization/acceleration  zone,  one  can  indirectly  measure  the  relative 
change  in  the  length  of  the  ionization/acceleration  zone  as  a  function  of  background  pressure.  For  a  thruster  like  the 
HiVHAc  EDU2,  the  part  of  the  ionization/acceleration  zone  that  is  inside  the  discharge  channel  is  roughly  the  same 
as  the  erosion  zone.  With  a  reasonable  assumption  about  how  the  upstream  edge  of  the  ionization/acceleration  zone 
is  anchored,  the  breathing  mode  frequency  can  be  used  to  determine  how  the  upstream  edge  of  the  erosion  band  is 
shifting  (the  downstream  edge  is  at  the  exit  plane  of  the  thruster).  This  information  can  be  used  to  establish  an  upper 
limit  on  acceptable  background  pressure  for  which  the  change  in  erosion  band  length  is  considered  acceptable.  The 
same  information  can  also  be  used  to  quantify  the  change  in  erosion  band  length. 

Figure  16  shows  a  plot  of  the  normalized  length  scale  as  a  function  of  normalized  background  pressure  generated 
from  high-speed  camera  data.  Normalized  length  scale  is  the  ratio  of  calculated  length  scale  to  the  predicted  length 
scale  at  zero  background  pressure.  For  example,  a 
normalized  length  scale  of  0.8  means  that  the  calculated 
length  scale  is  80%  as  long  as  the  length  scale  at  zero 
background  pressure.  This  figure  shows  a  very  consistent 
decrease  in  length  scale  as  the  background  pressure  is 
raised. 

To  establish  background  pressure  limit  for  future  test,  a 
few  assumptions  will  be  made.  For  a  typical  Hall  thruster, 
the  acceleration/ionization  zone  length  is  a  few  tens  of 
millimeters.  ’  ’  Assuming  a  few  millimeters  change  in 

the  acceleration/ionization  zone  length  correspond  to  a  few 
millimeters  change  in  the  erosion  band  length  and  assuming 
the  design  margin  on  the  erosion  band  length  was  chosen  to 
be  a  few  millimeters,  one  can  establish  0.9  as  the  limit  for 
how  much  change  in  normalized  length  scale  is  considered 
acceptable.  Using  Fig.  16,  the  normalized  pressure 
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Figure  16.  Normalized  length  scale  calculated 
from  high-speed  video  versus  normalized 
background  pressure. 


corresponding  to  a  normalized  length  scale  of  0.9  is  about  2-3.  Normalized  pressure  of  2.5  translates  to  about  1.2e-6 
Torr/(mg/s).  In  other  words,  the  erosion  band  length  will  have  changed  by  at  most  a  few  millimeters  if  the  operating 
pressure  is  less  than  1.2e-6  Torr/(mg/s).  Note  that  this  criterion  does  not  account  for  changes  in  plasma  energy  and 
overall  erosion  rate,  only  erosion  band  length.  For  reference,  the  criterion  on  background  pressure  based  on  a 
maximum  of  2%  increase  in  thrust  due  to  ingestion  is  5e-7  Torr/(mg/s).^16 

B.  Updated  Ingestion  Model 

In  a  previous  work,  a  two-stream  ingestion  model  was 
proposed  to  account  for  ingestion  of  background  neutrals  in 
a  Hall  thruster.^  The  part  of  the  model  that  predicts  the 
properties  of  the  external  stream  turned  out  to  yield 
unphysical  results  leading  to  the  conclusion  that  linearly 
super-imposing  the  internal  (from  propellant  manifold)  and 
external  (from  facility)  streams  on  top  of  each  other  does 
not  accurately  model  the  ingestion  physics.  Nevertheless, 
calculating  the  ingestion  rate  using  the  change  in  discharge 
current  should  be  a  valid  approach  since  every  ingested 
neutral  will  produce  an  electron  upon  ionization.  The  ion 
will  be  accelerated  through  a  part  of  the  potential  difference 
between  the  anode  and  the  cathode,  and  the  electron  will  Figure  17.  Illustration  of  the  simplified  two- 
travel  the  other  part  of  the  potential  on  its  way  to  the  anode,  stream  ingestion  model.  I4  Ving  indicate 

Figure  17  shows  an  illustration  of  the  two-stream  ingestion  voltages,  not  distances. 
model.  For  the  present  study,  the  model  is  modified  so  that 

the  ingestion  rate  can  be  calculated  without  making  any  assumptions  about  the  ingestion  area.  In  the  previous  study, 
the  ingestion  area  was  assumed  to  be  the  same  for  all  background  pressures  at  each  throttle  point. ^ 

The  updated  two-stream  ingestion  model  still  assumes  two  streams  of  neutrals.  The  first  stream  comes  from  the 
propellant  manifold  while  the  second  stream  comes  from  the  facility.  The  streams  have  the  same  species 
composition,  propellant  utilization  efficiency,  and  current  utilization  efficiency.  Under  these  assumptions,  the  ratio 
of  discharge  current  to  mass  flow  rate  for  the  ingested  stream  can  be  assumed  to  be  equal  to  that  of  the  main 
propellant  stream.  This  ratio  can  be  calculated  using  performance  data  from  the  single-string  integration  test.  The 
discharge  current  is  then  plotted  against  the  background  pressure  and  a  linear  fit  is  performed.  Zero-background- 
pressure  performance  of  the  thruster  and  ingested  current  can  then  be  calculated.  Ingested  mass  flow  is  then 
calculated  from  the  ingested  current  using  the  aforementioned  discharge-current-to-mass-flow-rate  ratio.  Table  4 
shows  the  measured  thrust,  discharge  current,  and  calculated  ingestion  rate.  Note  that  one  entry  in  the  calculated 
ingestion  is  negative.  This  unphysical  data  point  is  likely  the  result  of  a  poor  curve-fit  or  measurement  uncertainty. 
Recall  the  accuracy  of  the  mass  flow  controller  is  ±1%  of  reading. 
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Electron  path 

Channel  j  !  I  j 
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_ Table  4.  Summary  of  ingestion  rate  calculation. _ 

Disch.  Disch.  Thrust,  mN  Discharge  current,  A  Ingestion  rate  in  percentage 

voltage,  power,  _ of  supplied  mass  flow  rate 


V 

kW 

lx 

3x 

5x 

lOx 

lx 

3x 

5x 

lOx 

lx 

3x 

5x 

lOx 

300 

1.5 

85 

90 

- 

95 

4.95 

5.06 

- 

5.31 

0.5% 

2.4% 

- 

6.9% 

300 

3.0 

188 

191 

- 

196 

10.03 

10.06 

- 

10.32 

0.5% 

0.8% 

- 

3.3% 

400 

2.0 

105 

109 

- 

no 

5.04 

5.26 

- 

5.24 

-0.9% 

3.1% 

- 

2.8% 

400 

3.2 

172 

175 

- 

181 

8.02 

8.05 

- 

8.19 

0.3% 

0.6% 

- 

2.3% 

500 

2.0 

92 

95 

96 

- 

4.06 

4.18 

4.19 

- 

0.4% 

3.1% 

3.3% 

- 

500 

3.0 

139 

144 

146 

149 

5.99 

6.08 

6.13 

6.22 

0.0% 

1.5% 

2.3% 

3.7% 

500 

3.9 

182 

186 

188 

193 

7.74 

7.77 

7.90 

7.96 

0.3% 

0.7% 

2.3% 

3.1% 

C.  Dimensional  Analysis  on  Rising  Pressure  and  Shortening  lonization/Acceleration  Zone 

While  data  from  both  the  high-speed  camera  and  the  high-speed  Langmuir  probe  suggests  the 
ionization/acceleration  zone  is  shortening  and  the  downstream  boundary  is  receding  into  the  channel  with  rising 
background  pressure,  the  root  cause  is  unclear.  We  now  perform  a  dimensional  analysis  in  an  attempt  to  find  the  root 
cause. 
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Equation  (2)  shows  the  relationship  between  the  plasma  potential  and  charge  density  in  the  form  of  Poisson’s 
equation.  In  this  equation,  O  is  the  plasma  potential,  p  is  the  charge  density,  and  Sq  is  the  permittivity  of  free  space. 
The  corresponding  scaling  equation  is  shown  in  Eq.  (3).  Note  that  the  Laplacian  of  O  scales  as  V/L^  where  V  is 
voltage  and  L  is  length  scale. 


Li  ~ 


V^®  =  -p/eo 
1  1 

-^MRa  +MRi 


(2) 

(3) 


In  Eq.  (3),  Li  is  the  characteristic  length  scale,  is  the 
electron  density,  MRa  is  the  anode  mass  flow  rate,  and  MRi 
is  the  ingested  mass  flow  rate.  Eigure  18  shows  a  plot  of  the 
normalized  length  scale  predicted  by  the  above  dimensional 
scaling  against  ingestion  percentage  for  the  data  set  in  the 
present  study.  This  figure  shows  that  the  results  under¬ 
predict  the  change  in  length  scale  by  about  one  order  of 
magnitude.  A  more  complex  approach  is  needed  to 
accurately  predict  the  change  in  transport  that  occurs  with 
rising  background  pressure. 

VII.  Conclusions 


0%  1%  2%  3%  4%  5%  6%  7% 


The  present  study  combined  high-speed  camera  and  Ingestion  Percentage 

high-speed  Langmuir  probe  data  to  present  a  novel  picture  Predictions  of  normalized  length  scale 

of  the  effect  of  facility  background  pressure  on  the  plasma  ^  function  of  ingestion  percentage, 
oscillations  in  the  HiVHAc  EDU2  Hall  thruster.  In  the  process  of  studying  these  oscillations,  we  discovered  that  the 
discharge  of  the  HiVHAc  EDU2  at  the  operated  conditions  is  dominated  by  breathing  mode  oscillations.  Spokes  can 
occasionally  form  but  are  quickly  extinguished. 

The  high-speed  camera  data  showed  that  the  breathing  mode  frequency  varies  from  10  to  36  kHz  and  tends  to 
rise  with  rising  background  pressure.  At  the  same  time  the  magnitude  of  the  light  intensity  oscillation  inside  the 
discharge  channel  falls  with  rising  background  pressure.  Data  from  the  high-speed  Langmuir  probe  showed  that  the 
fluctuation  in  plasma  density  and  temperature  rises  with  background  pressure  10  cm  downstream  of  the  exit  plane. 
Data  from  the  two  diagnostics  combined  suggested  that  the  ionization/acceleration  zone  of  the  thruster  was 
shortening  and  receding  into  the  discharge  channel  with  rising  background  pressure.  If  true,  the  decreased 
characteristic  plasma  length  would  correspond  to  a  reduction  in  the  amount  of  field  lines  that  electrons  need  to  be 
transported  across  inside  the  discharge  channel  and  may  lead  to  a  corresponding  decrease  in  transport  being 
provided  by  the  breathing  mode.  This  phenomenon  would  explain  the  decrease  in  the  magnitude  of  light  intensity 
oscillations.  At  the  same  time,  the  increase  in  distance  between  the  downstream  boundary  of  the  acceleration  zone 
and  the  cathode  created  an  increased  reliance  on  the  breathing  mode  to  carry  the  electron  current  from  the  cathode  to 
the  channel,  which  gave  rise  to  an  increase  in  fluctuation  magnitude  in  plasma  properties  as  seen  by  the  Langmuir 
probes. 

The  scaling  relationship  between  breathing  mode  frequency  and  characteristic  length  scale  was  then  used  to 
establish  a  maximum  acceptable  background  pressure  for  which  the  erosion  band  length  will  have  changed  by  at 
most  a  few  millimeters.  This  maximum  acceptable  pressure  was  found  to  be  L2e-6  Torr  per  mg/s  of  anode  mass 
flow  rate. 

An  update  to  the  simple  two-stream  ingestion  model  was  made  to  decrease  its  dependency  on  unknown 
quantities.  Although  the  model  is  still  not  capable  of  accurately  predicting  the  properties  of  the  external  stream,  it 
can  be  used  to  calculate  the  rough  ingestion  rate  without  constraining  the  effective  ingestion  area. 

A  simple  dimensional  analysis  was  performed  regarding  how  the  ingestion  of  background  neutrals  can  cause  the 
ionization/acceleration  zone  to  shorten  and  recede  into  the  discharge  channel.  Although  the  resulting  prediction  is 
qualitatively  in  the  right  direction,  the  prediction  is  off  by  roughly  one  order  of  magnitude.  Eurther  work  is  needed  to 
determine  the  underlying  mechanism  that  links  rising  background  pressure  with  shortening  ionization/acceleration 
zone. 

A  number  of  assumptions  were  made  in  the  course  of  analyzing  the  data  presented  in  this  paper.  The  most 
crucial  assumption  was  that  the  scaling  equation  for  the  breathing  mode  frequency  is  correct.  An  important  next  step 
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will  be  to  determine  the  validity  of  said  scaling  relationship.  To  experimentally  verify  the  breathing  mode  frequency 
scaling  relationship,  one  can  deploy  three  diagnostics  simultaneously.  The  first  diagnostic  will  be  the  high-speed 
camera,  which  will  be  used  to  confirm  the  presence  of  the  breathing  mode  and  measure  the  associated  frequency. 
The  second  diagnostic  is  laser-induced  fluorescence  velocimetry,  which  will  be  used  determine  the  length  and 
location  of  the  acceleration/ionization  zone  as  well  as  the  velocity  of  the  ions.  The  third  diagnostic  is  either  a  set  of 
thermocouples  or  a  calibrated  thermal  camera  to  extract  the  anode  temperature,  which  will  be  used  to  calculate  the 
neutral  velocity.  The  background  pressure  will  then  be  varied  at  different  discharge  voltages  and  currents  to  produce 
the  data  needed  to  study  the  breathing  mode  frequency  scaling  relationship. 
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